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A new semi-analytical approach for numerical solving of Cauchy 
problem for functional differential equations 

Josef Rebenda * Zdenek Smarda t Yasir Khan * 


Abstract 

One of the major challenges of contemporary mathematics is numerical solving of various problems 
for functional differential equations (FDE), in particular Cauchy problem for delayed and neutral 
differential equations. Recently large variety of methods to handle this task appeared. In the paper, we 
present new semi-analytical approach for FDE’s consisting in combination of the method of steps and a 
technique called differential transformation method (DTM). This approach reduces the original Cauchy 
problem for delayed or neutral differential equation to Cauchy problem for ordinary differential 
equation for which DTM is convenient and efficient method. Moreover, there is no need of any symbolic 
calculations or initial approximation guesstimates in contrast to methods like homotopy analysis 
method, homotopy perturbation method, variational iteration method or Adomian decomposition 
method. The efficiency of the proposed method is shown on certain classes of FDE’s with multiple 
constant delays including FDE of neutral type. We also compare it to the current approach of using 
DTM and the Adomian decomposition method where Cauchy problem is not well posed. 


Introduction 


For the purpose of clarity, we consider the following functional differential equation of n-th order with 
multiple constant delays 

u {n \t) = . ,it ( " _1) (£), u i(f- ri),u 2 (t- r 2 ),.. ,,u r {t-T r )), (1) 

where ujf — Tj) = (u(t — Ti ), u'(t — Tj),..., — t,)) is mj-dimensional vector function, to, < n , 

r 

i = 1, 2,..., r, r G N and /: [£ 0) oo) x R n x is a continuous function, where u = m i- 

»=l 

Let t* = max{ri, T 2 ,..., T r }, m = max{mi, TO 2 ,..., m r }, m < n. In case m = n equation (JT]) is of neutral 

type, otherwise it is delayed differential equation. Initial function needs to be assigned for equation (Q]) 

on the interval [to — t*,to]. Furthermore, for the sake of simplicity, we assume that <f>(t) G C n ([to — £*Ro])- 

Investigation of equation m is important since there is plenty of applications of such equations in real life. 

As examples, we mention models describing behaviour of the central nervous system in a learning process, 
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species populations struggling for a common food, dynamics of an autogenerator with delay and 
second-order filter, systems controlled by PI or PID regulators, evolution of population of one species etc. 
For further models and details, see e.g. jT|. 

In the last two decades, various methods such as homotopy analysis method (HAM) [2], !3j, homotopy 
perturbation method (HPM) [4], variational iteration method (VIM) [Sj, Adomian decomposition method 
(ADM) 0, H, i, Taylor polynomial method [5J, Taylor collocation method uni and differential 
transformation method (DTM) 11 . jT2], [13] . [ Til have been considered to approximate solutions of 
certain classes of equation m in a series form. However, in several papers, for example Evans and 
Raslan [6], Mohammed and Fadhel [14] initial problems are not properly defined. The authors use only 
initial conditions in certain points, not the initial function on the whole interval, thus the way to obtain 
solutions of illustrative examples is not correct. Moreover, transformation formulas used in the calculations 
are very complicated. 

The crucial idea of our concept is to combine differential transformation method and general method of 
steps (more details on method of steps can be found for instance in monographies Kolmanovskii and 
Myshkis [T or Bellen and Zennaro ng). This approach enables us to replace the terms involving delay 
with initial function and its derivatives. Consequently, the original Cauchy problem for delayed or neutral 
differential equation is reduced to Cauchy problem for ordinary differential equation. Also the ambiguities 
mentioned above are removed. Further, while ADM, HAM, HPM and VIM require initial approximation 
guess and symbolic computation of necessary derivatives and, in general, n-dimensional integrals in 
iterative schemes, presented method is different: Cauchy problem for FDE is reduced to a system of 
recurrence algebraic relations. 

1 Preliminaries 


This section contains several well-known results which we compare to our approach later. 

The concept of differential transformation was used by Zhou |16] who applied it to solve linear and 
nonlinear initial value problems in electrical circuit analysis. 

Differential transformation of the k —th derivative of function u(t) is defined as 


U(k) = 


r d k u{t ) 


k\ 


dt k 


( 2 ) 


t=t 0 


where u(t) is the original function and U(k) is the transformed function. Inverse differential transformation 
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of U(k) is defined as follows: 


OO 


u{t) = Y,U{k){t-t 0 ) k . ( 3 ) 

k=o 

Many transformation formulas were derived from definitions ([2]) and Jd]) in the past. We mention several 
formulas for illustration: 

Theorem 1 ( m, v) and vi) from definition) Assume that F(k), G(k), H(k) and Ui(k), 
i = 1 ,..., n, are differential transformations of functions f(t), g(t), h(t ) and Ui(t), i = 1 ,. .., n, 
respectively. Then 


i) 

If 

^ g 

II 

then 

F(k)= { -±±^G(k + n). 

ii) 

If 

f(t) =g(t)h(t), 

then 

m=EtoG(l)H(k-l). 

in) 

If 

m = t n , 

then 

F(k) = 6(k — n), 6 is the Kronecker delta symbol, to 

iv) 

If 

< 

II 

-to 

then 

F(k) = 


v) 

If 

f{t) = sint, 

then 

F(k) = 

f(-l)V^ ifk = 2n+ l,n £ No t _ 0 

10 if k = 2 n, n £ No 

vi) 

If 

f{t) = cos t, 

then 

F(k) = 

k = 2n,n eN 0 ^ _ Q 

[^0 if k = 2n + 1, n £ Nq ’ ° 

vii) 

If 

/w = nr=i«i(*). 

then 





k k—s i 

k—s i — • 

■■Sn-1 




m=J2 E-- 

E Ui(si)... U n -\{s n -x)U„{k Si ••• s n ). 


Sl=0s2=0 Srx—0 


Theorem 2 (mi, m, mi) Assume that F{k), G(k) are differential transformations of functions fit), 
gift), where a > 0 is a real constant. If 

f{t)=g{t-a), then F{k) = ^(-l) J_fe j a l ~ k G(i), N -*• oo. (4) 

Using Theorem [2] and the formula i) in Theorem [l] differential transformation formula for function 
f(t) = 4nrg(t ~ a ) can be deduced. 

Theorem 3 ( [12] ) Assume that F(k), G(k) are differential transformations of functions f(t), g(t), 
a > 0. If 

m = ^g(t-a), then F(k) = £ (-1)*"*"" L * ) N -+ oo. (5) 

i=k+n ^ 
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Using Theorems [2J [3] and formula vii) in Theorem [Q any differential transformation of a product of 
functions with delayed arguments and derivatives of that functions can be proved. However, such formulas 
are complicated and not easy applicable for solving functional differential equations with multiple constant 
delays (see for example m , Ha, m- 

2 Main results 

Consider equation (ID subject to initials conditions 

u{t 0 ) = uo,u\t 0 ) =ui ,... ,u (rl-1) (fo) = u n -1 (6) 

and subject to initial function (f>(t ) on interval [to — t*,to] such that 

<t>(to) =u(to),(j>'(to) =u'{t 0 ),...,^ n - 1 '>(t 0 ) =u ( ' n ~ 1) {t 0 ). (7) 

First we apply the method of steps. We substitute the initial function and its derivatives in all places 
where unknown functions with deviating arguments and derivatives of that functions appear. Then 
equation m changes to ordinary differential equation 

u (n) {t) = f(t, it(i), u'(t), u (n_1) (f), it - ri), $ 2 (t - r 2 ),..., & r {t ~ T r )), (8) 

where &i(t — n) = ( <j>(t — n), — Tj),..., ^ mi \t — r,)), rrii < n, i = 1,2,..., r. Now applying DTM we 

get recurrence equation 

{k+ k ^ [ U{k + n) = r(k, 17(0), [/(l), ...,U(k + n-l)y (9) 

Using transformed initial conditions and then inverse transformation rule, we obtain approximate solution 
of equation (ID in the form of infinite Taylor series 

OO 

u{t) = Y^umt-tof 

fe =0 

on the interval [to, to + a], where a = min{ri, r 2 ,..., r r }, and u(t) = (j>{t) on the interval [to — t*, to]- We 
demonstrate potentiality of our approach on several examples. 

Example 1. Consider the problem that was solved using Taylor series method by Sezer and 
Akyuz-Dascioglu J9j and using DTM by Arikoglu and Ozkol [12| . 

u"{t) - tv!(t - 1) + u(t - 2) = -t 2 - 2t + 5, (10) 
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u( 0) = —1, u'{— 1) 


- 2 , 


—2 < f < 0. 


( 11 ) 


First, we remark that such formulation of problem is not correct since if we take for instance t = —1, then 
u(t — 2) = u(— 3) is not defined at all. If we omit condition —2 < t < 0 in (fill) , then according to Sezer and 
Akyuz-Dascioglu [9] we are looking for solutions of a problem involving mixed conditions. However, this is 
not a Cauchy problem, thus it is not clear what kind of solution are we looking for since uniqueness of 
solution is not guaranteed. Let us see consequence of this fact. 

Applying current approach of using DTM on both sides of equation (fTUl) and conditions (fill) . Arikoglu and 
Ozkol [T2] obtained recurrence relation 

(k + l){k + 2)U{k + 2)~Y, E (fc-fcr + 1)( 1 )(-l) hl -* +fcl - 1 ^(/i 1 )5(fc 1 -l) 

fci=0 hx=k-ki+l ' 1 ' 

N ,, s 

+ E ( fc ) {~2) hl ~ k U{hi) = - S(k - 2) - 2S(k - 1) + 5 S(k) (12) 

and transformed mixed conditions 


N 

u(o) = —i, Y.kum-i?- 1 = -2. 

fc=0 

Taking N = 4 and using inverse differential transformation formula (Jd)), they claimed solution 
u(t) = — 1 + t 2 which has the same form for all TV > 4. It can be easily verified that this u(t) is exact 
solution of m, cd for all f eR. Nevertheless, this coincidence is possible only for solutions in expected 
polynomial form. Generally, since TV is a finite number, we only get approximate solutions. 

Applying different approach using Taylor series method, Sezer and Akyuz-Dascioglu [9] obtained 
one-parameter class of solutions in the form u(t) = t 2 — 1 + a\{t + 2), where ai 6 1 is a parameter. It is 
easy to check that, again, such u(t) is exact solution of (flUl) . (fill) for all t € R. Furthermore, this class of 
solutions contains the solution achieved by Arikoglu and Ozkol jT2l . 

Now recall that to formulate Cauchy problem for (flUl) correctly, we have to prescribe initial function <p(t) 
on [—2,0] satisfying conditions © as well: 


u(t) = <j>(t) for t £ [—2, 0), u(0) = </>(0) = uq, u'( 0) = ^(0) = u±. (13) 


Thus a solution of Cauchy problem HDD, CD then should be expected in the form 


u(t) 


ip(t), t £ [0,1] 

m , [-2,0] 
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It is obvious that neither current approach of using DTM nor the other approach is suitable for solving 
Cauchy problem cud, m since, in general, there are infinitely many initial functions satisfying conditions 
nm> (see m below as an example), whereas those approaches always give the same result. 

On the other hand, our approach enables to find unique solution of Cauchy problem (fTUll . (THU) . For initial 
function cf>(t) = — 1 + t 2 satisfying (1111) we have simple recurrence relation 


U(k + 2) = 


2 *( fc ) fc>0 
(fc + l)(fc + 2) ’ 


From this relation and the fact that U( 0) = —1, 17(1) = 0, we get 17(2) = 1, U(k) = 0 for k > 3. Hence the 
exact solution is u[t) = — 1 + t 2 , generally on [0, oo), and this solution of Cauchy problem (fTOl) . (fldl) is 
unique. 

Further, for initial function 

<t>{t) = -t 2 - At - 1 (14) 


that also satisfies conditions CD, using DTM combined with the method of steps we obtain different 
solution of m, ci- In this case we have recurrence relation 


c (t+ 2 ) = k>o. 


(15) 


(fc + l)(fc + 2) 

From initial conditions u(0) = 0, u'(O) = —4, we have C/(0) = 0, 1/(1) = —4 and from (fl5ll we get 17(2) = 1 
17(3) = —2/3, 17(4) = —1/6, U(k) = 0 for k > 5. Thus we have exact solution in the form 


u(t) 


—1 — 4 t + t 2 — 1 1 3 — ^f 4 , t e [o, l] 
—t 2 — 4t — l, te [-2,0] 


and, again, this solution is unique. However, it is not possible to obtain this solution using either current 
DTM approach or the other mentioned method. 

Example 2. Consider Cauchy problem consisting of first order differential equation with one constant 
delay 

u'(t ) = — it(f) — -u(t — a) + a, a > 0, (16) 

a a 


and initial function 

cA(t) =t 2 , t £ [—a, 0]. 


(17) 


Obviously, u(t) = t 2 is unique solution of Cauchy problem (fTTH) . (1171) on [0,oo). From initial function (fTTl) 
we deduce initial condition «(0) = 0, which is transformed by differential transformation to 17(0) = 0. 
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First, we apply current approach of using DTM represented by Theorem [2) Transformed equation ( 1161 ) has 


the form 

1 1 N 

ck +1 )u(k +1) = -u(k) - - y (-i y~ k 
a a tt 

If we try to solve (fl8l) for IV = 1, we get 


17(1) = 17(1) + a 


a l ~ k U(l) + aS(k). 


(18) 


(19) 


for k = 0. It implies that there would be a solution only for a = 0 which is in contradiction to our 
assumption that a > 0. 

Next, if we solve m for TV = 2, we have 

17(1) = -17(1) - - [17(0) - aU (1) + a 2 17(2)1 + aS(0), (20) 

a a 

for k = 0, which after some rearrangements gives 


0 = — aU (2) + a, 


( 21 ) 


hence 17(2) = 1. Calculation for k = 1 does not provide any new information, while for k = 2 we get 
17(3) = 0. It can be easily computed that for TV > 2, we always get 17(2) = 1 and U(k) = 0 for k > 3. 

Note that we do not have any information about 17(1) for TV > 2. We can interpret it such that 17(1) may 
be an arbitrary constant C £ R. Taking this into account, we deduce that in this case solution u(t) is in 
the form 

u{t) = t 2 + Ct, C g K. (22) 

It is not difficult to verify that all such functions indeed are solutions of m satisfying initial condition 
u(0) = 0. It means that using current DTM approach we obtained one-parameter family of solutions of 
equation (ITU with initial condition u( 0) = 0, not a unique solution of Cauchy problem (fl6l) . (1171) . since we 
did not utilize initial function at all. Furthermore, from m we can see that the result of applying current 
DTM approach also may depend on the choice of TV. 

On the other hand, applying the new approach we derive relation 

(fc + \)U(k + 1) = -U(k) - - \S(k - 2) - 2 aS(k - 1) + a 2 6(k) 1 + aS(k), (23) 

a a 

which can be simplified to 

(k + 1)17(A: + 1) = ^17(fc) - i, 5{k - 2) + 2 S(k - 1). (24) 
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Transformed initial condition is the same as in the previous case, 17(0) = 0. For k = 0,1,2,... we have 


U( 1 ) = -U{ 0 ) 

2U (2) = it/( l) + 2 

3C7(3) = -C/(2) — - 
a a 


^[/(l) = 0, 
— 1!7 ( 2 ) = 1 , 


->E/(3) = 0, 


which gives unique solution 


u(t) = t 2 , t G [0, a]. 


(25) 


Applying the classical method of steps, we can see that (1751) is unique solution corresponding to Cauchy 
problem (ITU . (1T71) on [0, a], thus the new approach of using DTM is in perfect agreement with 
well-established results. 

Example 3. Consider delayed differential equation of the third order 


u'"{t) = —u(t) — u(t — 0.3) + e t+0 ' 3 


(26) 


subject to the initial function 


and conditions 


< f>(t ) = e 4 , t < 0 


u( 0) = 1, 
«'(0) = -i, 
u"(0) = 1. 


(27) 


(28) 


This problem was solved using the Adomian decomposition method (ADM) by Evans and Raslan [6] and 
later using current DTM approach by Karakoc and Bereketoglu jT3j and again using ADM by 
Blanco-Cocom et al. |7J. 

Straightforward observation gives the information that, as Blanco-Cocom et al. [7j point out, it is enough 
to consider only (1^51) and (1271) . since conditions (1751) are not independent of initial function tp(t) defined in 
m- However, in fact, in all mentioned papers authors did not use initial function (1271) at all, they solved 
problem (1751) . (1751) which is not a Cauchy problem. Karakoc and Bereketoglu [13] tried to rectify the 
situation of not using (1271) by excluding this condition from formulation of the studied problem. 


Unfortunately, this step led to the same curiosity observed in Example 1, when for instance u(— 0.3) is not 
defined. In any of the cases, uniqueness of solution is not guaranteed. 

In both papers using ADM the authors obtained approximate solution using iterative scheme containing a 
triple integral and compared the result to function u(t) = e _t which is a solution of (l26l) . (l28l) and satisfies 
(EH as well. Karakoc and Bereketoglu m solved equation (l26l) using current DTM approach without the 
dependence on the initial function and determined recurrence relation 

N fh \ 1 

(k + 1 )(k + 2)0 k + 3)U(k + 3) = -U(k) - ^ (-l) h '~ k ( 1 j (0 ,3) hl ~ k U(h 1 ) + -(-l) fc e°' 3 , (29) 

h\=k ^ ' 

The authors solved (l29l) for TV = 6,8,10 and compared obtained approximate solutions to solution 
u(t) = e _t . 

In contrast to complicated formulas mentioned above, our approach gives simple recurrence relation 


U(k + 3) = 


-U(k) 


From initial conditions 


(k + l)(fc + 2)(k + 3) ’ 
and recurrence relation (1301) we have 


k > 0. 


(30) 


17(0) = 1, 17(1) = -1, 17(2) = i, £7(3) = -±, 17(4) = i ..., U(k) = 


Using inverse differential transformation ([3| we obtain a solution of (l26l) . (l27l) in the form 

00 (— ~\\ k 

»« = E V- (t =< 31 ) 

k -0 

It is the closed form unique solution of Cauchy problem (l26l) . (l27l) which cannot be reached using either 
ADM or current approach of using DTM as described in above mentioned papers and EE only 

approximation of the solution is possible. 

Now consider initial function 

0(t) = ^t 2 -t+ 1, t€ [-0.3,0]. (32) 

This function satisfies conditions (l28l) as well. However, applying classical method of steps shows that 
function 


u{t) 


^ 0.3 


2.23 


e"* + 




0.3 i 

+ -- 1 2 + 1.37 - 1.39 

o z 


/—0.3 



(33) 
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is unique solution of Cauchy problem (1^1) . (15^1) on interval [0,0.3], nevertheless it is completely different 
from (ED). Using of our approach leads to recurrence relation 


(k + 1 )(fc + 2 )(jfc + 3 )U(k + 3) = -U(k) - ls(k - 2 ) + 1.3 6{k - 1) - 1.39 5(k) + e 0 3 (34) 

2 k\ 

which together with conditions (1281) transformed to 17(0) = 1, 17(1) = —1, 17(2) = 7 enables to calculate 


_9 89 -I- e 0 3 9 3 — e 03 —1+e 0 3 

£7(3) = ——— 57 ——> 17(4) = —— 77 -—j £7(5) = 1 + 6 


£7(6) = 

£7(9) = 


3! 

2.39- 2 e 0 ' 3 


6 ! 

-2.39 +3 e 1 
9!~ 


, £7(7) = 


4! 

-2.3 + 2 e 0 - 3 
7! 


£7(8) = 


5! : 

1 - 2 e 0 3 
8 ! : 


(35) 


0.3 


Precise investigation of this sequence reveals that U(k), k = 0,1, 2,... are coefficients of Taylor expansion 
of function (l33l) in neighbourhood of 0, thus the calculation resulted in an approximation of unique solution 
of Cauchy problem (|26l) . (l32l) on [0, 0.3] as expected. Note that in this case, not even approximation of 
solution is possible using either ADM or current DTM approach as described in papers 0 , m and m- 


Example 4. The following Cauchy problem for first order neutral differential equation 

u'(t) + — ?/(i — 1) = u(t ) + u{t — 1), t> 0, (36) 

u(t) = —t, forte [—1,0]. (37) 

was investigated by Fabiano m- The author used semidiscrete approximation scheme to approximate 
unique solution of problem m, ED- The exact solution which can be calculated by method of steps is 
given by 

u(t) = t- i + ^e*, t e [0, i]- (38) 

We remark that the sewing condition 

^(O) + ^ , (-l)=^(O) + 0(-l) 

is not fulfilled hence the derivative of solution of problem E3, ED is not continuous at 0. 

Applying DTM combined with method of steps, we obtain relation 

{k + 1 )U{k + 1) - ^S(k) = U(k) - S(k - 1) + S(k), (39) 
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which together with initial condition u(0) = 0 acquired from (1371) and transformed to 17(0) = 0 implies 


U(1) -\ = U(0) +1 
217(2) = 17(1) — 1 
317(3) = 17(2) 
417(4) = [7(3) 


= 

-+U( 2) = 
^[/(3) = 
^[/(4) = 


5 

4’ 

1 _ 1 

8 - 4-2!’ 

1 _ 1 

24 “ 4-3!’ 
1 

4 ^ 4 !’ 


If we write 17(0) and 17(1) as 17(0) = ^ | and 171) = + 1 and perform inverse transformation, the 

solution can be expressed as 


u{t) = 


4-0! 


1 +iV+-Att7 2 + 


4-1! 


4-2! 


1 , 1 4 

t + ——r t + • ■ 


4-3! 


4-4! 


4 • k\ 


t K + ••• = t-- 


1 “ 1 


E ^ ( 4 °) 


4 ^ k\ 

k—0 


which is Taylor expansion of exact solution (1381) . Hence using our approach we are able to identify unique 
solution of Cauchy problem (1361) . (1371) in closed form, which is not the case in Fabiano’s paper El- 


Conclusion 


• We conclude that combination of method of steps and differential transformation method (DTM) 
presented in this paper is powerful and efficient semi-analytical technique suitable for numerical 
approximation of a solution of Cauchy problem for wide class of functional differential equations, in 
particular delayed and neutral differential equations. No discretization, linearization or perturbation 
is required. 

• There is no need for calculating multiple integrals or derivatives and less computational work is 
demanded compared to other popular methods (Adomian decomposition method, variational 
iteration method, homotopy perturbation method, homotopy analysis method). 

• Using presented approach, we are able not only to obtain approximate solution, but even there is a 
possibility to identify unique solution of Cauchy problem in closed form. 

• A specific advantage of this technique over any purely numerical method is that it offers a smooth, 
functional form of the solution over a time step. 

• Another advantage is that using our approach we avoided ambiguities, incorrect formulations and 
ill-posed problems that occur in recent papers, as we observed in examples. 
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• Finally, a subject of further investigation is to develop the presented technique for equation © with 
state dependent or time dependent delays. 
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